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The Gross-Pitaevskii (GP) equation is a long-wavelength approach widely used to describe the 
dilute Bose-Einstein condensates (BEC). However, in many physical situations, such as higher den- 
sities, this approximation unlikely suffices hence one might need models which would account for 
long-range correlations and multi-body interactions. We show that the Bose liquid described by the 
logarithmic wave equation has a number of drastic differences from the GP one. It possesses the 
self-sustainability property: while the free GP condensate tends to spill all over the available volume 
the logarithmic one tends to form a Gaussian-type droplet - even in the absence of an external trap- 
ping potential. The quasi-particle modes of the logarithmic BEC are shown to acquire a finite size 
despite the bare particles being assumed point-like, i.e., the spatial extent emerges here as a result 
of quantum many-body correlations. Finally, we study the elementary excitations and demonstrate 
that the background density changes the topological structure of their momentum space which, in 
turn, affects their dispersion relations. Depending on the density the latter can be of the massive 
relativistic, massless relativistic, tachyonic and quaternionic type. 
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1. INTRODUCTION 

The Gross-Pitaevskii (GP) approximation is a long- 
wavelength theory commonly used for describing dilute 
Bose-Einstein condensates (BEC), such as trapped al- 
kali gases. As a matter of fact, the GP equation (also 
known as the cubic Schrodinger equation) is the result 
of few physical assumptions [l|: the mean- field approx- 
imation, neglected excited states and multi-body inter- 
actions (three and more particles scatter at one point), 
two-body interaction being assumed of the contact type 
(delta- function) , neglected anomalous contributions to 
self-energy when using the perturbation theory [2j . There 
is no warranty whatsoever that these assumptions must 
hold for all physical examples of BEC. Indeed, higher 
densities (common for the Bose liquids made of noble 
gases or nuclear matter under certain conditions Q), low- 
dimensional effects or the presence of mixtures can under- 
mine their validity and the theory needs drastic modifica- 
tions of a universal form. Such modifications are usually 
made in the form of corrections to the GP equations or, 
alternatively, to the expression for the energy functional: 
while the energy in the GP approximation is quartic in 
|\&|, *S? being the wave-function of condensate, the correc- 
tions can be of the sextic and octic power, some examples 
being given in Ref. jfj|. Yet, even those corrections may 
not suffice, e.g., if higher densities lead to the long-range 
correlations and multi-body interactions. In that case 
one should account for all the powers of |^| which im- 
plies the usage of the non-polynomial functions which will 
necessarily appear in the wave equations for condensate. 



An example of the usage of non-polynomial functions for 
certain class of condensates has been demonstrated, for 
instance, in Ref. Q, although the known examples refer 
to effectively lower-dimensional models so far. 

In this paper we heuristically introduce the conden- 
sate described by the non-linear Schrodinger equation 
of a non-polynomial kind, namely, the logarithmic one 
(LogSE): 



T-L - /3 _1 ln(a 3 |*| 2 ) * = 0, 



(1) 



where ^ refers in general to the complex-valued wave 
functional, /3 and a are the parameters of the theory, 
and % is the operator the form of which is determined 
by a physical setup - for instance, in a non-relativistic 
theory % — H — ihd t where H — J^- + V cxt (x) is the 
Hamiltonian operator. The physical motivation behind 
this equation as well as its unique properties are listed in 
the Appendix. Since the pioneering works Q this equa- 
tion received much attention - its applications were found 
not only in the extensions of quantum mechanics, but 
also in quantum optics (3), nuclear physics Q, transport 
and diffusion phenomena [9( , open quantum systems and 
information theory (Io| . effective quantum gravity and 
physical vacuum models [TlL fT3| . Thus, a theory of Bose 
liquids can be yet another interesting area of application 
of LogSE, certainly worth being explored. 

To begin, let's suppose that we have the BEC of N 
identical particles whose wave-function is described by 
the equation (JT|) which in our case becomes 
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*(2,i)=0, 
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where m is the mass of the condensate particle, 14 x t(x, t) 
is the trapping potential. Here the notations for /3 are as 
follows: our j3 is equivalent to 1/6 from Ref. [6| and has 
a sign opposite to /3 from the corresponding part in Ref. 

The particle density of the condensate is defined as 
usual: 



n = \^(x,t)\ 2 



(3) 



where <j>(xi, t) describes the single-particle state of the zth 
boson in the fully condensed state. 

For now, we leave aside the question of the microscop- 



ical structure which might lead to Eq. be it long- 
range correlations, multi-body scattering, fermionic de- 
grees of freedom or anomalous energy contributions, and 
study the logarithmic Bose liquid as it is (although, some 
elements of the microscopic theory are given below in the 
section allotted to the Bogoliubov excitations) . On top of 
the properties which directly follow from those listed in 
the Appendix, this type of condensate has certain unique 
properties as a supcrfluid. 

First of all, let's consider the effective potential energy 
for the logarithmic BEC. If one defines the Lagrangian 
density as 
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Vextl*! 2 -/?- 1 !*! 2 [ln(a 3 |<5| 2 )-l] 



(4) 



then 
ft (7 



the 

dL 



corresponding 
+ V- 



Eulcr-Lagrange equation, 
= 0, yields Eq. ©. 
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Thus, the effective potential energy density is given by 



Va = -p- x n [ln(na 3 ) - l] 



(5) 



and thus it opens up (down) and has the local non-zero 



minima (maxima) at n c 



for the negative (posi- 



tive) p. In what follows we call this potential logarithmic 

- due to the property dVp /dn cx In (na 3 ) which yields the 
logarithmic term in the Schrodinger equation. 

This potential is non-analytic yet regular at n = 

- while the logarithm itself diverges there, the factor n 
recovers the regularity, see Fig. [TJ In fact, the poten- 
tial always has the Mexican-hat shape if plotted as a 
function of "J. This leads to the drastically new fea- 
tures which will be discussed below. It is not diffi- 
cult to check also that the Ginzburg-Landau potential 
(the quartic potential supplemented with the chemical- 
potential quadratic term) is one of the perturbative lim- 
its of the logarithmic one - if, for instance, one ex- 
pands Vp in the vicinity of its local extrema n c ^t- Vp — > 
-/3" 1 [(a 3 /2)n 2 -n- l/(2a 3 )] + C[(na 3 - l) 3 ]. Obvi- 
ously, this expansion only makes sense when the density 
is close to the special value 1/a 3 but in general the po- 
tential (J5J) can not be replaced by the polynomial one and 
thus requires a non-perturbative treatment. 

For instance, as long as the potential (0 vanishes in 
the extrema n xt, they describe the equilibrium config- 
urations for "free" particles. These equilibria can be 
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FIG. 1: The field-theoretical potential ([5]) of the logarithmic 
BEC in the dimensionless units: —/3a 3 (vertical axis) and 



-3/2 



(horizontal axis). 



stable or unstable and here appears something which 
might look counterintuitive at first sight. The negative 
P means that the the field-theoretical potential opens up 
and the above-mentioned extrema are minima but the en- 
ergy functional is not bounded from below Q , therefore, 
the stable equilibrium solution can only exist in a finite 
volume or in the additional external potential ("trap") 
[l3| . On the contrary, when /3 is positive then the poten- 
tial ([5]) opens down but the energy functional is bounded 
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from below. An apparent contradiction can be resolved if 
one recalls that we are dealing with the normalized wave- 
function and thus its square (the condensate density n) 
cannot grow arbitrarily large. The physical consequence 
is that the coherence length is always larger than the ef- 
fective size of the droplet. In Sec. [2] and [3] we will show 
that this leads to one of the most important features of 
the Bose liquid with logarithmic nonlinearity: at posi- 
tive (3 it can be self-bound and stable (as to form the 
droplet-like object) even without any trapping potential. 

Second, the logarithmic term in Eq. ([2]) can be associ- 
ated with the certain kind of quantum-informational en- 
tropy oc — J | \T/ 1 2 In (Sl^l 2 )^ 3 ^, measuring the degree 
of quantum spreading of the BEC as a collective quan- 
tum object, as being described in the Appendix. The 
only difference from the interpretation presented there 
is that in the BEC case the logarithmic term describes 
the quantum information transfer between the states of 
the correlated particles which form the many-body state 
whereas the temperature is a formal measure of the 
collisionless interactions inside the Bose liquid at zero 
(thermal) temperature, see the Appendix. 

Third, using the Madelung representation of the wave 
function and the hydrodynamic form of the wave equa- 
tion ([5} [3, EEj], the zero-temperature (collisionless) 
equation of state of the logarithmic BEC in the leading 
order with respect to the Planck constant is described by 
a kind of the Clapeyron-Mendeleev law, 



p — Po = — (m(3) n + 0(h 2 ) oc T^n. 



(6) 



For comparison, the equation of state for the CP conden- 
sate would be p oc n 2 , thus, the logarithmic Bose liquid 
is more "ideal" than the Gross-Pitaevskii one yet non- 
trivial. From here it follows that the logarithmic BEC is 
the only one where the speed of sound (phonons) does not 
depend on the density in the leading order with respect 
to the Planck constant: 



l/y/m\j3\+0(H). 



(7) 



Meanwhile, the speed of sound for the GP condensate 
scales as a square root of density. Due to this property 
the logarithmic BEC can be used to describe or mimic 
the induced relativity and several high-energy and grav- 
itation al p henomena, along the lines described in Refs. 
[IH [3 [13 > although considering these topics would bring 
us too far from the scope of the current paper. 

Finally, in the GP case it is possible to differentiate 
between the "attractive" and "repulsive" interaction de- 
pending on a sign of the coupling constant. In our case 
the logarithm changes sign when density goes across the 
value a~ 3 which results in changing the sign of the inter- 
action in Eq. ([2]) and the type of the condensate depends 
not only on the sign of (3 but also on whether the value 
of na 3 is smaller or larger than one. 



2. THE LOGARITHMIC BEC IN THE EMPTY 
SPACE 

While many of the properties discussed in this section 
hold also for the condensate in a trap, it is more clarifying 
to see them on the example of the released condensate. 
If it were the GP-type Bose gas then its properties would 
be nearly trivial: once the trap is off the GP conden- 
sate spreads uniformly all over the available space (if it 
is repulsive), or shrinks down until the new phase arises, 
in the case of attraction. The behavior of the free log- 
arithmic BEC will be shown to be much more complex, 
mostly due to the properties of the logarithm mentioned 
above. 

If we assume the isotropic case for simplicity and work 
in the center-of-mass reference frame of the condensate 
then ty(x,t) = &(\x\) exp (— ifj,t/K), fi being the chemical 
potential, and the equation ([5} reduces to the ordinary 
differential equation for the modulus of wave function 
\& = v/ n(r) as the function of the radial coordinate only, 



2m 



A* + [n + /T 1 In (a 3 * 2 )] * = 0, 



(8) 



where A*? = r~ 2 d r (r 2 d r ^) . Further it will be conve- 
nient to introduce a value of dimensionality of the length 
squared, Ap = h 2 /3/2m which measures the strength of 
the logarithmic interaction in terms of length scales, and 
a value of dimensionality of the energy, E a = H 2 /(2ma 2 ) 
which is the de Broglie energy corresponding to the 
length scale a. Upon applying the normalization con- 

oo 

dition, J nd 3 x = A-k J ^> 2 r 2 dr = A, we deduce that the 
o 

normalizable solution is only possible for positive j3. The 
exact expressions for, respectively, the particle density 
and the chemical potential are: 



n(r) = * 2 =n e- {r/a ^\ 



(9) 



^ = 3(3- 



1 - In 



A 1 / 3 a 



-p- x \n{N/N ), (10) 



where n Q = A/(7r 3 / 2 a 3 3 ) = (e/a) 3 A/A , a p = \f~A~p, = 
h\J j3/2m, and Ao = (e^/irap / a) 3 is the number of par- 
ticles at which the chemical potential changes its sign. 
The energy functional for the logarithmic BEC can be 
defined in a standard way (see, e.g., Sec 6.1 of the book 

My- 

oo ^ 

£[n]=47r J (V^n) 2 -/?- 1 [In (na 3 ) - l] j nr 2 dr, 

o 

(11) 

which yields the expression for the energy per particle: 



E/N = 3/T 1 



In 



A 1 / 3 a 



/r^i-MA/Ao)] 

(12) 
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One can see that £ as a function of N vanishes in the 
origin and at e-/V , and has a positive- valued local maxi- 
mum, i.e., it is bounded from above (see Fig. [2J 



E/N < 0- 



(13) 



where the value of N at which the maximum occurs 
is given by the transcendental equation In (N/N ) + 
Nq/N = 1. Notice that for the non- negative values of 
E there exists some sort of degeneracy - one value of 
energy corresponds to two values of N. This will be dis- 
cussed below in more detail. Further, if we treat E as a 
function of the coupling then it obviously vanishes when 
P^ 1 approaches zero. However, it can also vanish at the 
nonzero value of the coupling which can be made obvious 
by rewriting the expression as 



E/N = --p- 1 \n {p-'/Po 1 ), 



(14) 



where Pq 1 = Tre 8 / 3 N~ 2 / 3 E a represents the non-zero 
value of the parameter P at which this energy vanishes. 
From this expression one can also see that the energy is 
bounded from above 



E(P) < ^e- 1 ^ 1 = ^ne' 



' 3 N^E a , (15) 



with this local maximum occurring at P m ^ — P /e — 

ne 5/3 N -V3E a . 

It is important to note that both the chemical poten- 
tial ([TO]) and the ground-state energy (|T2"T) depend on the 
parameter a as lna^/a and as we show below this term 
can only be positive as ap > a. 
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FIG. 2: The total energy (solid line) and the chemical poten- 
tial (dotted line) for the logarithmic BEC with positive /3 in 
absence of any trap as a function of the number of particles. 



2.1. Effective potential, stability and critical mass 

For the analysis of the solutions of the generic equation 
of the form 



-^v 2 + ^ cxt (f, i) + F(|*| 2 )- M 



* = 



it is often convenient to introduce the notion of the ef- 
fective potential V Q g(r) = V cxt {x, t) + F(\^>(r)\ 2 ) which 
in our case is defined as 



V eff (r) = -/T 1 ln(a 3 <I> 2 (r)). 



(16) 



Being evaluated on the solution (|9]) it becomes the 
quadratic function of radius- vector: 



V eS (r) 



1 



-muj eS r +fi 



-fkU, 



off, 



where 



Wcff : 



2p- 1 /h. 



(17) 



(18) 



Therefore, the effect from the logarithmic term here is 
that it effectively introduces the trap with the frequency 
which is controlled by P whereas the system can be 
viewed as a particle of energy /i whose center of mass 
is located in this trap. Then the stability of the solution 
against a small perturbation is ensured by the fact that 
this frequency is non-negative if P is non-negative. 

This is also confirmed by the Vakhitov-Kolokolov sta- 
bility criterion [l9] which in our case is equivalent to the 
non-positivity of the derivative of the chemical potential 
with respect to N. From Eq. (fTOf one can derive that this 
derivative equals to ~P~ 1 /N and thus it is non-positive 
indeed. 

There is, however, an additional issue of stability. 
From Eqs. (| 12[) or (|14[) it is clear that for the non- 
negative values of E there exists some sort of degener- 
acy - one value of energy is achieved at two values of N 
or P, respectively. This indicates that the system spon- 
taneously 'chooses" either of these values at no cost of 
binding energy, and the so-called mechanical instability 
occurs {20]. Thus, our condensate is impeccably stable 
only at 



E < 0, 



(19) 



which means that the number N must be bounded from 
below 



N > N n 



eNn. 



(20) 



The physical meaning of the latter is that one cannot 
create the stable logarithmic condensate with just any 
amount of initial (bare) particles - it must be larger than 
the value iV m i n . In other words, in the presence of certain 
long-range correlations and multi-body scattering the ini- 
tial Bose system must have the mass larger than the crit- 
ical one, mN m i n , for the Bose condensation to happen. 
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2.2. Self-sustainability and droplet formation 

An important thing to notice is that indeed the den- 
sity of the free condensate with the positive (3 does not 
become the uniform one - in empty space the conden- 
sate rather tends to take the shape of a spherical Gaus- 
sian droplet, with a characteristic size of order ap. Our 
droplet, however, is different from the ordinary ones. As 
long as its density n(r) is described by the Gaussian the 
major portion of mass and energy is localized in a finite 
volume - similarly to the liquid droplets in the clas- 
sical world. The difference is that the Gaussian vanishes 
at the infinity only, therefore, our droplet is essentially 
inhomogeneous and has no sharply defined external sur- 
face, hence, its stability is ensured by the nonlinear quan- 
tum effects in the bulk rather than by the surface tension. 

Further, from Eqs. (HUJ) and (fH]) it follows that for the 
Gaussian solution the following relation holds 

V-E/N = -r X , (21) 

which essentially means that the chemical potential 
counted from the energy per particle is negative in our 
case. A similar feature was observed in real superfluids, 
such as 4 He [2l|, therefore, models involving the logarith- 
mic nonlinearity might find an application there. 



N do not contribute to AE and can be taken either to 
#(°) or to S4 0) while 5* oc N\nN: 

— Tvj/iSiJ;, (23) 

S* = k B NlnN, (24) 
E^ =N/3- 1 (1 + In N ), (25) 

where = (ksP)^ 1 is the associated temperature. 
Thus, we can formally divide the energy into two types 
- into the entropy and the part which is linear w.r.t. N. 
Therefore, the latter is linear w.r.t. density too, while 
J5(°) /N can be considered as the energy of a free quasi- 
particle. Alternatively, the iV-linear part can be consid- 
ered as a reference entropy = N(l + In Nq) which is 
supposed to be positive. 

In general, such definition of entropy can be deduced 
for any shape of the logarithmic BEC (e.g., in the trap) 
as the integral in Eq. (|2"2"j) can be always presented in the 
form N In N + Nf where / is an N- independent function 
of the parameters of a system. If we require the positive- 
ness of the additive part of the energy (it means that / is 
always positive) the entropy part of the energy Eq. (|2"5)) 
determines its lower limit at given N. 

2.4. Emergent spatial extent 



2.3. Quantum Shannon entropy and non-thermal 
temperature 

Let us consider now the notion of the information en- 
tropy introduced in the Appendix. In the context of the 
Bose-condensation phenomenon the entropy is a mea- 
sure of information stored in the microscopical configura- 
tion of N bosons which form the droplet where the square 
of wave function determines the density value and con- 
stant 3 defines the information entropy reference frame. 
Besides, becomes a thermodynamic-like measure of 
the quantum spreading of the droplet as a whole. 

To start with Sq, , it is natural to measure the constant 
S in the units of a 3 : S = Z]J N a 3 , where Zq is a dimen- 
sionless constant. Then the entropy (IA.2I) can be written 
as 

k^Sy = k^S^ + J n\n(a 3 n)d 3 x, (22) 

where = ks^nZo is the reference entropy. One 

can physically justify the identification of an entropy 
using the fact that the energy functional is bounded 
from below. The reasoning is the following. The en- 
ergy needed to cut the A^-particle droplet into two equal 
parts is AE = — N(3~ 1 ln2, into three parts it is AE = 
—Nj5~ 1 In 3 and so on. Therefore, the maximum binding 
energy which is saved in a droplet and can be realized is 
no more than AQ = In N. So we can use the ther- 

modynamic definition of entropy through AQ — T$ AS 1 *. 
The parts of the energy (TT21 which are proportional to 



From Eq. © one can see that while the logarithmic 
condensate is "made" of N bare particles which are as- 
sumed to be point-like the resulting Gaussian droplet ac- 
quires a finite size ap. Besides, it is clear that the par- 
ticle density of the logarithmic condensate is bounded 
from above, n ^ no , which implies that the average vol- 
ume occupied by a condensate's quasi-particle is bounded 
from below 

(Vol) > {^a fi f/N. (26) 

Together with the arguments coming from the analysis of 
elementary excitations, see the paragraph after Eq. ([55]), 
this means that the effective size of the quasi-particle can 
be non-zero at the non-zero P~ l ■ Recalling the meaning 
of ap as an effective size of the Gaussian droplet, we see 
that the latter can accommodate no more than N max _ 
quasi- particles, i.e., 

N < A> max w ivra^/(Vol) = N/(6y/Z), (27) 

where tilded values refer to the condensate's quasi- 
particle modes. 

However, the inequalities of this type are trivial in a 
sense; they follow entirely from the normalization con- 
dition and thus they are not necessary specific of the 
logarithmic model. In our case, the better bound can 
be derived from physical arguments. It is clear that the 
iV-linear part (|25|) must be non-negative because this en- 
ergy can only increase when N does. Thus, we have to 
impose IniVo ^ — 1 which is equivalent to riga 3 ^ e 4 N or 
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ap a/(7r 3 / 2 e 4 / 3 ) (we remind that ap is the mean size 
of the Gaussian droplet). As long as n ^ we obtain 

na 3 < n Q a 3 Ne 4 , (28) 

which implies that a r /r ^ 1, where a r = a/(A r e 4 ) 1 / 3 is 
the quasi-particle effective size and r is the mean inter- 
particle distance. Thus, the smallest volume of the sys- 
tem of N particles is nonzero but of the order a 3 . This 
is crucially different from what the GP approach would 
yield - there no such effect appears, the Gaussian-type 
droplet can only form when the gas is placed into the 
external trapping potential of a special kind. It is inter- 
esting that the spatial extent exists in the worldvolume- 
formulated theories [22[ (where it is postulated from the 
beginning), and also it is known to arise as the quan- 
tum effect in the quantum mechanics of point particles 
on non-commutative spaces [23j . In our case, however, no 
spatial extent was initially postulated (the bare particles 
were assumed to be point-like, according to the standard 
quantum-mechanical approach) and no quantum commu- 
tation rules were modified - instead, the non-zero size is 
an emergent quantum phenomenon due to the peculiar 
type of correlations between particles' wave-functions. 

2.5. Effective potential and finite volume 

As shown above, the value a 3 is the smallest possible 
volume of the system of N particles whereas a r ~ a/N 1 ^ 3 
is the size of a quasi-particle (for brevity, in this section 
we neglect the 0(1) factors such as those containing e and 
7r). One can argue that the condition (|28|) means that 
states with na 3 > N can not be physically reached as 
particles commence to "overlap" each other. Therefore, 
the effective potential, 

V(n) = Vp- [in, (29) 

is confined to the given bounds, see Fig. [3] In par- 
ticular, it is naturally bounded from below because 
na 3 In (na 3 ) <C N\nN. A similar condition holds for 
the central density which implies that ap/a > 1 (the 
droplet 's volume must not be smaller than the smallest 
one). Therefore, our potential for the positive j3 has two 
minima, one being located at zero, and the other being at 
the limit value corresponding to maximally dense pack- 
ing. The latter is usually deeper and hence defines the 
ground state. For the system with a finite number of par- 
ticles, such as our droplet, the density of a dense packing 
is much larger than its central density no. The depth of 
this minimum can be controlled with a quadratic term, 
see Fig. [3J If one considers this potential with the ex- 
tra quadratic (chemical potential) term one can find that 
negative chemical potentials shift the local maximum to 
the large values. At ln(na 3 ) = 1 — [i m axlfi~ X the min- 
imum at zero becomes deeper and the maximally dense 
packing state is not the ground state anymore hence the 
system can not be self-sustained anymore. 
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FIG. 3: The effective potential <[29j) for the logarithmic BEC 
for zero and large negative chemical potentials 

The size of the droplet is defined by a parameter /3 and 
by putting more and more particles into the system the 
maximum density can be reached when na 3 « N which 
means ap a. One might expect that at a given N one 
can get a denser condensate by putting it inside a trap, 
but the size of this trap should be smaller than ap. In 
this case the zero minimum of the effective potential is 
responsible for the ground state. It should be noted here 
that if one selects the harmonic trap to identify such a 
phase transition it might be rather difficult as the ground 
state for both phases is still described by the Gaussian 
functions, only of a different width. A clearer way would 
be to use a non-harmonic trap so that a phase transition 
should reveal itself through the change in the shape of 
the droplet. 

As for the case of the negative (3, the ground state is 
determined there by the minimum of the effective poten- 
tial at na 3 = 1. We found that this state can not be 
self-bound and can be confined only in a trap. 



3. LOGARITHMIC BEC IN A TRAP 

First we consider the logarithmic BEC placed in the 
infinite-walled spherical square-well potential of radius 
d. As one can see from ([31]) , the condensate density in a 
trap can be much larger than for a free droplet, see Fig. 
[4] for the case of positive f3. When the size of trap d is 
much larger than the characteristic size of the droplet ap, 
the condensate does not "sense" the trap and behaves 
like a self-bounded Gaussian droplet. As long as d ~ 
ap the condensate begins loosing its Gaussian shape and 
becomes denser in the center. In the case of the negative 
(3 there is no localized solution and condensate tends to 
fill the whole trap of any size. It should be noted that the 
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condensate density profile does not depend on the quasi- 
particle's size and only determines the possible number 
of particles in the condensate. 
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FIG. 4: The density of the logarithmic BEC in the infinite- 
walled spherical well of radius d. 

Now, let us place the logarithmic condensate into the 



J 



isotropic harmonic trap of frequency uj. Then the equa- 
tion ([2|) reduces to the following one 



2m 



A* + [|mw 2 r 2 



Ai-/3 _1 ln(a 3 * 2 )] * = 0, (30) 



where a w = y^h/muj is the characteristic width of the 
harmonic trap. Unlike the previous (free) case, the ex- 
act normalized solution of this equation exists both for 
positive and negative /3's and still has the Gaussian form: 



n(r) = n Q (tu)e-^/ a ^ 2 



(31) 



where 



n {u) = Ti V^N/aju, ap u 
and 



2A r - 



1 - X ^/me^ 



(32) 



X = - sign (/?) y/l + (2Ap/al) 2 = - sign {fi) y/l + {hwpf 

(33) 

so that the values e u p — (1 — X)/ ft and a are never neg- 
ative. The normalization condition yields the allowed 
value of the chemical potential: 



1 - X - In 



ma 2 N 2 / 3 1 - x 



= -/?- 1 ln(iV/iVoH), 



(34) 



r 



where 



(35) 



and in the limit to — > at the positive /3 the expres- 
sions expectedly reduce to those from the previous sec- 
tion, N (uj = 0) = N . 

By analogy with the spherical- wall case one can show 



J 



that when fko ^> /3 (which is equivalent to a w < ap) 
the central density in the trap can be considerably larger 
than the central density of the droplet. But of course it 
can not be infinitely large because the condensate density 
has an upper bound 



Further, the energy of the system is given by 



E[n] = 4tt 
hence the energy per particle becomes 



2mn 



i^Vn) 2 + ^muj 2 r 2 - /T 1 [in (na 3 ) - l] j nr 2 dr, 



5 , f ma 2 N 2 ^ l-x 
Z- X - ln \— 2 — 



= p- 1 [l-ln(JV/iVb(w))] 



(36) 



(37) 



r 



The relation //^ — E^/N — —/3 1 holds, similarly to the case of the condensate without a trap, which means that 
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the chemical potential counted from the condensate en- 
ergy per particle is negative (positive) for the positive 
(negative) j3. The plot of the energy in /3~ 1 Nq units for 
the logarithmic BEC in an isotropic harmonic trap as 
functions of a number of particles in Nq units is given in 

Fig. G3 
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FIG. 5: The energy of the logarithmic BEC in an isotropic 
harmonic trap versus the number of particles for different val- 
ues of /hiv. 



3.1. Effective potential and linear stability 

The effective potential in this case is defined as 
1 



V e $(r) = -mu 2 r 2 - In (a 3 # 2 (r))- 

Being evaluated on the solution (|3~Tj) it becomes: 

1 3 
Vcs(r) = -muj 2 s r 2 + fa- -ttiJ e g, 

where 



(38) 



(39) 



Weff = W fi - r 1 (l + sign 08) y/l + (fi^/3) 2 ) /h. 

(40) 

Therefore, the effect from the logarithmic term here is 
that it replaces the trap frequency with the modified one 
which depends on /3 and it is always non-negative regard- 
less of the signs of ui or /3. 1 Thus, the coefficient at r 2 



1 This effect is similar to what happens in quantum mechanics on 
non-commutative spaces where the motion of the particle placed 
in a harmonic trap stays oscillatory but it is determined not by 
the bare frequency but by the modified one containing the non- 
commutativity parameter [23ll . 



is non- negative regardless of the sign of j3 either, hence, 
the Gaussian solution pip is stable both at a positive 
and negative j3. Moreover, it stays stable at arbitrary 
small Jvuj/3 unless the latter vanishes exactly - then the 
negative-/? solution decays, in agreement with the result 
from the previous section. This feature shows another 
difference of the logarithmic BEC from the GP one - the 
latter decays when the non-linear coupling (proportional 
to the scattering length) becomes negative. These prop- 
erties are essentially rooted in the sign-changing property 
of the logarithm discussed in the introductory part of the 
paper. 



3.2. Mechanical stability and related bounds 

As for the mechanical stability the arguments given for 
the no-trap case can be easily generalized for our case: 
for the logarithmic condensate to be stable in a harmonic 
trap the following inequality must hold 



^- 1 [l-ln(JV/JV o (w))]<0, 



(41) 



whose physical implications are to be discussed now. 

As long as the coupling is not necessarily positive now 
we must distinguish two cases: 

(a) Positive j3. The inequality (RTj) brings, as in the no- 
trap case, the lower bound for the number of condensate 
particles, with the only difference that now it becomes 
w-dependent 



2vr 



\X\ 



1 



a 

3/2 , , 

py e 2(3b <-■•>. 



(42) 



so it can be regulated by adjusting the trap's frequency. 
Similarly to the no-trap case, this bound implies that one 
can not create the logarithmic condensate with just any 
amount of initial (bare) particles - it must be larger than 
the critical value which is frequency-dependent 

now. If the frequency energy is much smaller than the 
logarithmic coupling, hui -C then the critical value 
is a quadratic function of the frequency: 



iV (+) = eN 



+ 0((M3) 4 ), (43) 



where No was defined in the no-trap case. If the fre- 
quency energy is much larger than the logarithmic cou- 
pling, Hlu 3> /3~ 1 , then the critical value grows exponen- 
tially with an increasing frequency: 



N (+) 



7r 3/2 (<Wa) 3 exp 1 



(44) 



which means that the trapping potential dominates over 
the logarithmic term in the wave equation, and the con- 
densate goes into another phase. 
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(b) Negative (3, There appears the upper bound for the 
number of condensate particles which is also w-dependent 



N < AT(-) = eN (u) = tt 3 / 2 (^) 3 e^ 3 ^ 



2tt 



W-i 



3/2 



e -2( 3 lx|-5)^ 



(45) 



This critical value can be also regulated by the trap's 
frequency. If the frequency energy is much smaller than 
the logarithmic coupling, fruj <C /3" 1 , then the critical 
value scales as the inverse cube of the frequency: 



(46) 



where 



„ = a*/y/\Ap\ = uj x y is another length 

scale. When the frequency vanishes, diverges which 
indicates that at negative /3 there is no Gaussian so- 
lution for the logarithmic condensate in empty space, 
as was mentioned in the previous section. If the fre- 
quency energy is much larger than the logarithmic cou- 
pling, fouj 3> /9 j then the critical value vanishes expo- 
nentially with increasing frequency: 



N (-) 



3/2 



(a u /a) 3 exp | 



(47) 



as the trapping potential dumps the logarithmic-term ef- 
fects and thus destroys the logarithmic condensate. 



4. EXCITATIONS IN THE LOGARITHMIC BEC 

Small fluctuations in quantum Bose liquids can be con- 
ditionally cast into two classes - collective modes and el- 
ementary excitations. The former describe the motion of 
the liquid or a large part thereof as a whole, the latter 
are localized and can be viewed as particle-like objects. 
Although, due to the particle- wave duality there is no ex- 
act border between these classes: for instance, a phonon 
which is a quantum of collective vibrational modes can 
be also described as an elementary excitation via the Bo- 
goliubov approach. 

The collective mode of the logarithmic BEC is gov- 
erned by the second-rank tensor ( "acoustic metric" ) cor- 
responding to the conformally-flat Lorentzian manifold 
hj. That result can be easily transferred to the con- 
densed matter physics provided that the fundamental 
velocity constant is assumed to be the speed of sound. 
Here we consider the other class, elementary excitations, 
for the logarithmic condensate in the empty space. To 
this effect we perform the second quantization by going 
from the wave-functions to the field operators via the 
Bogoliubov decomposition 



(48) 



where n is the density of the condensate (for simplic- 
ity we assume the latter being uniform which is a valid 



approximation at the length scales are smaller compared 
with the size of the droplet), 8^> is the operator of ex- 
citations, [6^f, S^ 1 ^] = 9, and throughout the paper the 
non-polynomial functions of operators are defined as the 
Taylor series of matrices, as usual. One of the simplest 
forms of the Hamiltonian can be easily deduced from Eq. 
(fl"j) assuming the normal-ordering approximation 



H 



dr 



2m 



/T 1 : * t *ln(a 3 * t *) 



(49) 

where the colon means the normal ordering of the oper- 
ators. Upon applying the Bogoliubov decomposition, we 
obtain 



H 



= [ £o - /r 1 

2 



^1 + In (a 3 n ))] 6&54> 

5& 2 ) , (50) 



+ <5* t2 ) 



when omitting non-quadratic terms; here we denote eo = 
p 2 /2m — ft, p being the absolute value of the momentum 
vector. Upon applying the Bogoliubov diagonalization 
procedure we obtain 

H — ftN = ePb, (51) 
e 2 = [e - /T 1 (1 + In (a 3 n ))] 2 - /T 2 , (52) 

where w and b are the creation and annihilation oper- 
ators of the quasi-particles with momenta p, [b, b<] — 
[S^f,Sw], defined as the linear combinations 



6* = 

<y*t 



Kb + Vk 2 



IP, 



lb + KbK 



where 



«2 



«1 



V2 VL z > 
eo-/?" 1 (l + ln (a 3 n )) , 



-1/2 



k 2 = /3-\ 

and we can write the dispersion relation for the excita- 
tions of the logarithmic BEC in the form 



e(p) 



^- - ft - /3- 1 (1 + In (<z 3 n )) 
Am 



/3~ 2 , (53) 



and also we assume /3 to be positive for the remainder 
of this section 2 , therefore, the value of the chemical po- 
tential can be borrowed from Eq. (|10[) . As compared 



2 The condensate with the negative /? in the empty space can not 
be described by the Gaussian, therefore, its treatment would 
be more complex and model-dependent. Although, one cannot 
exclude that some results of this section will also be applicable 
for the negative-/? sector of the model. 
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with the dispersion relation for the excitations in the GP 
condensate 



£gp(p) 



the dispersion for the logarithmic Bose liquid has a num- 



J 




ber of drastic differences. 

The first thing to notice is that the requirement of non- 
negative energy points out the existence of the forbidden 
region of momenta (where e becomes negative or even 
complex which indicates instability issues). The forbid- 
den momenta region is given by the inequality (we are 
neglecting the chemical potential for the moment) 



-1 In (a 3 no) < ( ^ /T 1 [in (a 3 n ) + 2] if ln(a 3 n )>0, 

V 2 ™ /forbid 

V 2 ™ /forbid 



(54) 
(55) 



r 



and this region only exists if the upper bound ln (a 3 no)+2 
is positive-valued. The latter condition is equivalent to 
ro/e 2 / 3 > a where tq = n is the average spacing be- 
tween quasi-particles. Therefore, the appearance of the 
gap in the momentum space is in fact another manifesta- 
tion of the emergent spatial extent phenomenon discussed 
above. 

Thus, depending on the value of the background con- 
densate density the momentum space of the model can 
contain forbidden regions and thus can have a non-trivial 
topology. To see the full picture, we introduce two aux- 
iliary values of dimensionality energy 

£l = /3- 1 ln(a 3 n ) + fx, (56) 
e 2 = e 1 +2/3- 1 >e 1 , (57) 



then the dispersion relation (|53j) can be rewritten in the 
form 



<~-2 



(58) 



which is more convenient when classifying dispersion rela- 
tions according to the topology of the momentum space. 
We distinguish the following cases: 

(i) £2 < 0. Then ei must be negative as well, therefore, 
energy vanishes nowhere and the topological structure 
of the momentum space is trivial. The quasiparticle's 
spectrum is fully gapped and equivalent to that of the 
massive relativistic particle moving in the 4D spacetime 
with the fundamental velocity constant determined by 
the speed of sound c: 



0{p 4 /m 2 



P 

2m* 



where 



A = 



2A 



ei + e 2 



2m 





— m 


ei + 


£2 




- £2 



<D(p 4 /m 2 ), 
(59) 



(60) 



To allow the excitations of this type, the background con- 
densate density must be below the first threshold value 



n < tit, 



where 



Tlx 



1 N 



N 



e 5 7r 3/2 a 3 



(61) 



(62) 



(ii) e 2 = 0. Then e 1 = -2p~ l and thus it is still 
negative, in the momentum space there appears a zero- 
energy point, at p = 0, so the topological structure is 
still trivial and this case can be treated as a limit of 
the previous one. The quasi-particle's spectrum in the 
vicinity of this point is equivalent to that of a phonon 
or a massless relativistic particle moving in a spacetime 
with the fundamental velocity constant Cq: 



2 2 2 

e =c Q p 



where the speed of sound is given by 

c = 1/ \/m/3 5 



(63) 



(64) 



so that Co = c s , cf. Eq. 0. To allow the excitations 
of this type, the background condensate density must be 
equal to the first threshold 



n 



tit- 



(65) 



(iii) C2 > 0, ei < 0. To allow the excitations of this 
type, the background condensate density must lie be- 
tween the first and second threshold values 



tit < no < e n-r, 



(66) 



then in the momentum space there appears a zero-energy 
sphere of radius p = ^/2me 2 - The positive-energy modes 
are only located outside it, the quasiparticle's spectrum 
in the vicinity of this sphere is equivalent to that of a 
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tachyonic particle moving in a spacetime with the funda- 
mental velocity constant c 2 : 

e 2 = p 2 c 2 2 + m\c\ + O ((p 2 - 2me 2 ) 2 ) , (67) 

where the speed of sound and imaginary mass are 

C2 = \l ' 2 n 61 , m 2 = 2i. — — — m. (68) 
V 2m V e 2 - £ i 

In particle physics such models are usually rejected as in- 
compatible with Lorentz symmetry but in our case this 
can not be a valid reason: this symmetry is not a fun- 
damental but rather an emergent phenomenon which ap- 
pears at the lowest levels of energy. 

(iv) ti > 0, ex = 0. This happens when the background 
density approaches the second threshold 

no = e 2 n,T, (69) 

then in the momentum space there appears the zero- 
energy sphere of radius p = \/2mt2 and the zero-energy 
point at p = 0. This leads to degeneracy and the appear- 
ance of two zero-energy states. The first mode is identical 
to the mode described in the previous item upon setting 
ei to zero. The other mode, the one with vanishing p, 
has the formal dispersion relation: 

e 2 = -^-p 2 +0( P *), (70) 

which is of the quaternionic or 4D Euclidean (50(4)- 
symmetric) type. However, this dispersion is only valid 
at p — 0, therefore, this mode is non-propagating. 

(v) e 2 > 0, ei > 0. This corresponds to the density 
growing above the second threshold 

n > e 2 n T , (71) 

then the momentum space gets split by two zero-energy 
spheres, of radii \J1mt\ and \/2~m€2, which leads to 
degeneracy and the appearance of two zero modes. 
Positive-energy modes form two bands, p < y / 2mei and 
p > \J2me2- The first zero mode, located near the inter- 
nal sphere y/2mei, is described by the dispersion: 

£ 2 = fl^2 p 2 + 2£i(£2 _ £l)+0 ( (p 2 _ 2m£i) 2) (?2) 

m 

so it belongs to the above-mentioned quaternionic type 
with the only difference that now its energy is positive 
and the momentum does not vanish. The other zero- 
energy mode is of the tachyonic type, it is equivalent to 
the one described by Eq. (jfJT)) . 

To conclude, the topology of the momentum space of 
the logarithmic BEC excitations depends on the back- 
ground value of the condensate's density. In its turn, 
depending on the topology of the momentum space their 
dispersion relations can cover all cases described by the 
symmetries of the special orthogonal groups of order 
four - massive relativistic, massless relativistic, tachy- 
onic and quaternionic. The tachyonic (quaternionic) 



mode arises when the momentum space contains the zero- 
energy sphere so that positive-energy states are located 
outside (inside) it. One can also notice that the dis- 
persion relations become more and more exotic as the 
background density grows. This is in fact another mani- 
festation of the above-mentioned emergent size effect. 



5. CONCLUSION 

In this paper we considered the toy model based on 
the logarithmic Schrodinger equation and applied it to 
describing the ground state and excitations of some hy- 
pothetical Bose condensate. That is why we did not spec- 
ify why and how such sort of nonlinearity might exist 
in real Bose liquids. We only pointed out that LogSE, 
despite its simplicity, reveals an interesting class of so- 
lutions and it can possibly be adapted for a description 
of strongly- interacting Bose liquids. The important fea- 
ture of the logarithmic BEC is that it can exist as a self- 
bound droplet since the logarithmic nonlinearity uniquely 
combines both an attractive and a repulsive interaction. 
But a self- sustainable solution may exist for the gener- 
alized GP equation in the case of a very large scattering 
length as well, where the combination of an attractive 
two-body interaction and a repulsive three-body interac- 
tion allows to manufacture self-bound quantum droplets 
24]. Though one may suppose that such systems with 
a large scattering length are not "safe" against TV-body 
interactions and, for example, a four-body interaction 
might destroy a droplet if it is too attractive while some 
high-order interaction might stabilize it. However, the 
logarithmic type of an interaction naturally has both the 
attractive and the repulsive parts which are formally nei- 
ther two-body nor any other ./V-body interactions. There- 
fore, LogBEC does not need to introduce any extra in- 
teraction terms in order to get a stable self-sustainable 
solution. 

We found that the self-bound ground state, Gaussian 
droplet, can only exist at a positive parameter j3 and its 
energy is bounded both from below and from above. It 
is provided by the fact that the potential energy can be 
both positive and negative whose behaviour is defined by 
its logarithmic dependence on density. The effective po- 
tential for such state is of the harmonic trap type whose 
frequency is defined by a parameter /3. The logarithmic 
BEC with a negative parameter /3 can not have the self- 
bound state and only exists in a trap. We found that in 
the harmonic trap the solution is always of the Gaussian 
type for the parameter /3 of any sign and value. 

We justified the existence of the characteristic quasi- 
particle's size which limits the maximum density and the 
ground state energy of the condensate. This character- 
istic size naturally arises if we directly require the posi- 
tiveness of the additive part of the ground state energy 
(El. 

Finally, we studied the elementary excitations and 
demonstrated that depending on the topological struc- 
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ture of the momentum space their dispersion relations 
can be of the massive relativistic, massless relativistic, 
tachyonic and quaternionic type. 
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Appendix: Logarithmic Schrodinger Equation 

There exist at least two ways of how the logarithmic 
Schrodinger equation (LogSE) can be introduced. The 
original one is based on the separability argument - the 
LogSE is the only local Schrodinger equation (apart from 
the conventional linear one) which preserves the separa- 
bility of the product states: the solution of the LogSE 
for a composite system is a product of the solutions for 
uncorrelated subsystems [6] . 

The second way is based on arguments which come 
from open quantum systems and quantum information 
theory [TcJ. It is relatively less known and thus de- 
serves to be mentioned here. Let us consider a multi- 
particle (sub) system whose dynamics is described by the 
Hamiltonian-type operator H. Besides, this subsystem 
is in contact with its environment so that there is an 
exchange of energy and information. The state of the 
system is described by the vector |\E'). If the Hamilto- 
nian does not depend on a wave function then in the 
Schrodinger coordinate representation we recover the lin- 
ear differential equation for vj/. 

However, in general the interactions between the parti- 
cles comprising the subsystem depend on the distribution 
\^>\ 2 of the particles in the configuration space. To deter- 
mine this distribution, i.e., to extract, transfer and store 
the information in a particular configuration of matter, 
one requires a certain amount of energy per bit, call it 
e. The information acquired upon measurement of the 
state is proportional to the logarithm of the probability 
of an outcome 'J, i.e., 



I* = -log 2 (S|tf| 2 ) = -ln(S|vI/| 2 )/ln2, 



(A.l) 



and the associated Shannon entropy of the subsystem is 
given by 3 



5* = -Cfcs(*|ln(S|*| 2 )|*), 



(A.2) 



3 As far as we know, this kind of entropy was introduced first 
in Ref. [25l | and subsequently rediscovered by several authors 
afterwards. 



where ks is the Boltzmann constant, Q is the sign func- 
tion chosen in such a way that stays positive. This 
entropy approaches minimum on delta-like distributions 
and maximizes on uniform ones, therefore, it can be used 
as a measure of "spreading" of the probability distribu- 
tion described by Here the normalization factor H 
defines a measurement reference for the entropy because 
for continuous systems the latter is not absolute. For in- 
stance, one could establish the reference entropy as that 
for a uniform distribution hence if the subsystem has a 
fixed volume and the states are box-normalized then 5 
equals this volume. 

It should be also mentioned that this entropy is closely 
related to the logarithmic Sobolev inequality (known by 
physicists as the Everett-Hirschmann uncertainty rela- 
tion) which was conjectured by Everett [251 ] . rediscov- 
ered by Hirschmann [26| and proven by Beckner [27j 
using the Babenko inequality [28|]. It has been shown 
that this uncertainty relation is stronger than the Heisen- 
berg one, among its other applications is the proof that 
the energy in logarithmic nonlinear quantum mechan- 
ics is bounded from below [(|. Essentially, the Everett- 
Hirschmann uncertainty means the following: suppose 
we have the probability densities in the d-dimensional 
space and its Fourier transform which are normalized to 
N, f\V(r)\ 2 d d r = J\^{k)\ 2 d d k = N, *(k) = J"*(x). 
Then the Everett-Hirschmann inequality reads 

*(r)| 2 In |*(r)| Vr - J |*(k)| 2 In |#(k)|Vk 

> d(l + \mr)N -2N\nN, (A.3) 

where the l.h.s terms are easily recognizable as S\& for the 
position and momentum space, respectively. To prove 
that this inequality is stronger than the Heisenberg un- 
certainty, one has to use Eq. (|A.3|) to show that 



J |*(r)| 2 ln|*(r)| 2 d d r 



— (r- (r)) 2 



- (l + 2/d)\nN\, (A.4) 



and similarly for \l/(k); here ( . ) means an average, as 
usual. 

Let us go back to the Hamiltonian. The above- 
mentioned energy e thus contributes to the Hamiltonian 
of the form 



H^H' = H- £ log 2 (S|v]/| 2 ), 



(A.5) 



and the formal (non-thermal) temperature which can 
be associated with this kind of entropy is given by 
T<s; = (kB/3)- 1 = (dE'/dS 9 ) s = e/(fc B ln2), where 
E' = (#JH'J£) = E + (T^Sy is the total energy of the 
system [lOl . [Llj . Rewriting e in terms of /3, we recover 
LogSE in our notations (fTJ). For stationary states one 
can write it in the form 



H - p- 1 ln(H|#| 2 ) # = 



(A.6) 
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whereas the free energy is given by E = ( 1 I f |H|\I f ) = 
E' — (T^Sty. Unlike the free energy, the energy T^S^, 
is engaged in handling the information 7$ and thus un- 
available to do dynamical work. 

The Schrodinger equations of this type are suitable for 
describing subsystems in which the information is not 
conserved but being exchanged with environment, and 
where one can introduce some sort of temperature and 
entropy. Therefore, they can not be naively applied to 
the systems without any kind of irreversibility hence the 
negative results of the experiments [2!| are not surpris- 
ing. On the other hand, some Bose liquids might be 
good candidates, as being studied in this paper. Another 
candidate would be the yet unknown theory of quantum 
gravity where the debates still continue about the black 
hole thermodynamics and a possible information loss [30| . 
To conclude, we write down the most important proper- 
ties of LogSE: 

• Separability of noninteracting subsystems (as in the 
linear theory): the solution of the LogSE for the 
composite system is a product of the solutions for 
the uncorrelated subsystems; 

• Energy is additive for noninteracting subsystems 
(as in the linear theory); 



• Planck relation holds as in the linear theory; 

• All symmetry properties of the many-body wave- 
functions with respect to permutations of the coor- 
dinates of identical particles are preserved in time, 
as in the linear theory; 

• Superposition principle is relaxed to the weak one: 
the sum of solutions with negligible overlap is also 
a solution; 

• Free-particle solutions, called gaussons, have the 
coherent-states form, and upon the Galilean boost 
they become the uniformly moving Gaussian wave 
packets modulated by the de Broglie plane waves; 

• Expressions for the probability density and current 
are the same as in the linear theory. 

All these properties except the last one and, perhaps, 
second last and third last ones, are unique to LogSE 
among all other local nonlinear Schrodinger equations. 
Besides, many of these features are pertinent to the lin- 
ear Schrodinger equation which makes the logarithmic 
one a "minimal" nonlinear modification in a sense. 
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